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We consider the effect of the Coulomb interaction in strained graphene using tight-binding approx- 
imation together with the Hartree-Fock interactions. The many-body energy dispersion relation, 
anisotropic Fermi velocity renormalization and charge compressibility in the presence of uniaxial 
strain are calculated. We show that the quasiparticle quantities are sensitive to homogenous strain 
and indeed, to its sign. The charge compressibility is enhanced by stretching and suppressed by 
compressing a graphene sheet. We find a reduction of Fermi velocity renormalization along the 
direction of graphene deformation, in good agreement with the recent experimental observation. 
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I. INTRODUCTION 

Graphene is a two-dimensional crystal of carbon atoms 
with enormous interest on its unique featureJ^. Graphene 
sheet has attracted considerable attention because of its 
unusual electronic properties^ which follow from chiral 
band states and because of potential applications'^. The 
low energy quasiparticle excitations energies in graphene 
are linearly dispersing, described by Dirac cones at the 
edges of the first Brillouin zone. Charge carriers move 
strictly in two-dimensions since the graphene layer is just 
one atomic monolayer thick. 

Graphene is amenable to ex tern al influences incorpo- 
rating mechanical def ormat iorP^. Its band structure 
does not change for realistic strains less than 159^P^. 
The influence of long-range strains on t he electronic 
properties is a unique feature of graphen4^. At low- 
energy spectrum, strains give rise to a pseudomagnetic 
field which is added to the momentum operators^'^ and 
thus a gauge field couples to electrons. The most evi- 
dence of the unusual way in which strains affect the elec- 
tronic states comes from scanning tunneling microscope 
measurements of the electronic local density of states of 
graphene grown on platinurrPi'. An average compression 
of 10% creates effective fields of the same order of mag- 
nitude with the value observed in experimentJ^^l. 

Tuning the dynamics of massless carriers by appro- 
priately designed strain patternJl^ opens the way for 
novel applications of graphen^^l Strain can be induced 
in graphene either intentionally or naturally. Uniaxial 
strain can be induced by bending the substrates on which 
graphene is elongated without slippage. Elastic responses 
have been measured by pushing a tip of atomic force 
microscopes on suspended graphene. The presence of 
ripples in graphene samples and their influence on the 
electronic properties are open problems in the field for 
instance graphene on top of Si02 or SiC surface experi- 
ences a moderate strain due to surface corrugations or 
lattice mismatch. 

The Fermi velocity is an essential quantity in graphene 
because all the observable quantities depend on it. 
For the Dirac electrons in pristine graphene, it was 



showrP^H2SI that interaction effects also become notice- 
able with decreasing density in that the velocity is 
enhanced rather than suppressed, and that the influ- 
ence of interactions on the compressibility and the spin- 
susceptibility changes sign. These qualitative differences 
are due to exchange interactions between electrons near 
the Fermi surface and electrons in the negative energy sea 
and to interband contributions to Dirac electrons from 
charge and spin fluctuations. Recent experiment J^SIII] 
have been able to measure the renormalized Fermi veloc- 
ity and they claimed that the Fermi velocity is no longer 
constant but increases by decreasing the electron den- 
sity. On the other hand, it was showrP21, from first princi- 
ples calculations, that the group velocities under uniaxial 
strain exhibit a strong anisotropy. As the uniaxial strain 
increases along a certain direction, the Fermi velocity 
parallel to it decreases quickly and vanishes eventually, 
whereas the Fermi velocity perpendicular to it increases 
by as much as 25%. It was also shown that the work 
function of strained graphene increases substantially as 
strain increases. 

The role of long-range electron-electron interactions 
in a uniaxial strain on undoped graphene, has been re- 
cently studiecP^ using the renormalization group the- 
ory. They showed that while for small interactions and 
anisotropy the system restores the conventional isotropic 
Dirac liquid behavior, at intermediate coupling both the 
anisotropy and interactions can flow, in the renormaliza- 
tion group sense, toward large values thus signaling the 
emergence of anisotropic excitonic states. On the other 
hand, the dependence of the electron velocity on peri- 
odic deformations of graphene has also been investigated 
by another groufP^. The authors showed that the Fermi 
velocity anisotropy corresponds to the anisotropy of the 
quasiparticle spectrum energy. 

The positions of the Dirac points can be shifted in 
the opposite directions by applying anisotropic strains. 
The reason for this is that the time-reversal symmetry is 
preserved by strain in graphend^. However, the transla- 
tional symmetry is broken by lattice deformation. 

Our aim in this work is to study the Coulomb inter- 
action effects in uniaxially strained graphene particulary 
the anisotropy of the Fermi velocity renormalization and 
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the charge compressibility within the exchange interac- 
tion. In the Hartree-Fock theory, the positions of the 
Dirac points are shifted in the opposite directions when 
the full Brillouin zone calculation is taken into account. 
Our theory for strain dependence of quasiparticle veloc- 
ity renormalization in interacting Dirac electron systems 
is motivated not only by fundamental many-body con- 
siderations, but also by application and potential future 
experiments in the strain engineering field. 

The paper is organized as follows. In Sec. II we 
introduce the formalism that will be used in calculat- 
ing strained ground state properties which includes the 
many-body effects by using Hartree-Fock ( HF) approxi- 
mation. In Sec. Ill we present our analytical and numer- 
ical results for the renormalized Fermi velocity and the 
charge compressibility in doped graphene sheets. Sec. IV 
contains discussions and conclusions. 
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II. METHOD AND THEORY 
A. Noninteractacting strained graphene 

The electronic bands in graphene arise from the hy- 
bridization of pz orbitals which are localized at each car- 
bon atom. A good description of the electron states is 
obtained by assuming that electrons can hop between 
nearest-neighbor atoms. We consider a two dimensional 
honeycomb lattice in which a Pz electron hops between 
lattice points and it is governed by the nearest-neighbor 
tight-binding Hamiltonian 




Ho 



Ujalbj+H.c . 

<ij> 



(1) 



where Ui and bj are fermion field operators in sublat- 
tices A and B, tij are hopping integrals between lattice 
points and will be different among different neighbors. 
The Fourier component of the Hamiltonian is written as 



^0 




(2) 



where the factor f{k) is — J2i U e > with being 
the hopping energy of two nearest neighbor atom. Notice 
that a lattice deformation changes the off-diagonal terms 
in the Hamiltonian which describe hopping between the 
two sublattices which make up the honeycomb lattice. 
The induced scalar potential created by strain is inten- 
sively suppressed by screening and thus we do not take 
it into account. 

The honeycomb lattice geometry, three nearest- 
neighbor vectors and the lattice primitive vectors of 
a pristine graphene are depicted in Fig. [l] The lat- 
tice vectors are d^f 



- (l,v^)v^ao/2 and 4°^ 



FIG. 1: (Color online) Top) Honeycomb lattice geometry 
incorporates three nearest-neighbor vectors and the lattice 
primitive vectors of pristine graphene. Bottom) The first Bril- 
louin zone of undeformed graphene incorporates two Dirac 
points with high symmetry. 



(-1, \/3)%/3ao/2 where gq = 1.42 A is the C - C equilib- 
rium distance. The nearest neighbor vectors are defined 
by 



%/3 _1 
2 ' 2' 
ao(0,l) 

V3 1. 

«o(-^,-2) 



(3) 



Reciprocal lattice basic vectors are also defined by 
df'^ ■ — 2TrSij . The Brillouin zone which contains two 
Dirac points is shown in Fig. [T] (bottom) in which the po- 
sition of point in this perfect lattice is K 



X. 



The nearest neighbor vectors in strained lattice are de- 
fined'^ by 5' i_i = Sf'^ -\- ^ij^^j and thus the modifica- 
tion of these distances distorts the reciprocal lattice and 
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it is easy to find the reciprocal vectors under strain and 
the results are 



b'l 
b'2 



27r 
3x/3d 



3(1 + eyy) - Vie^y, \/3(l + exx) - 'i^xy 
-3(1 + eyy) - V^exy, a/3(1 + exx) + S^xy 



(4) 



where d = 1 + Cxx + £yy + O(e^) and Eyx denotes a shear 
strain and being zero in the uniaxial strain, respectively. 
The tensor for uniaxial strain along x direction is 



1 

-I 



(5) 



where v — 0.165 is Poisson's ratio. Notice that we are 
interested in uniform planar tension cases. Most impor- 
tantly, the Dirac point which is a symmetry point, will 
be shifted to a new position. For homogeneous strain 
the position of Dirac point is defined by the condition 
f{K) — 0. After some straight forward calculations, the 
new position of the Dirac point is given by 
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Notice that the reciprocal vectors of unstrained 
graphene appear in Eq. ^ instead of the distorted vec- 
tors,'2£l and the reason of that originally comes from the 
tight-binding Hamiltonian and the definition of f{k). In 
order to find the new position of the Dirac point, we 
only need to evaluate the hopping integrals. In fact, the 
change of bond-length leads to different hopping inte- 
grals among neighbors. Since the carbon atoms are out 
of equilibrium distance, we setP^ 
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FIG. 2: (Color online) a): Contour plot of the energy dis- 
persion for strained and interacting graphene sheet and b) 
a zoom of the energy dispersion around A'o-point. Notice 
that the center of box indicates the position of K'^ ( shift 
of the Dirac point in strained graphene without the electron- 
electron interaction) and clearly Kn moves to the left. These 
results are calculated by considering e — 0.05, aee = 0.5 and 
n =10^^cm~^. The deviation from the linear dispersion rela- 
tion associated to the trigonal warping effect is clearly shown 
in (a) using the full expression of f{k). 



9.1 ( l^'il 1-1 , , 

U = toe '0 ( ^0 1) , (8) 

where = 3.09ey gives rise to vp = lO^m/s. This 
definition satisfies a conditiorPSl in which dti/da = —6.4 
eVA^^. The dispersion relation is thus given by i?o(fc) — 
±|/(fc)|. For the sake of completeness, we thus expand 
the form factor around the /iT^-point to obtain its low- 
energy expression and it results in 

i i 

(9) 

From now on, we drop (0) index of the nearest-neighbors 
vectors for simplicity. Equation (|9]) is thus simplified as 



|/(/4 + q}\' « 4{qAf + tliql^f + tl{q.5^f 

+ [qA){q.5:i){tl-tl~tl) 

+ {qA){ql2){tl~tl-tl) 

+ {q.52){q.h){tl-tl^tl) (10) 

Therefore, the band-structure with arbitrary hopping in- 
tegrals is given by 

El = \f{KD +q)\^ = h^vlql + vlql) + c\xqy , (11) 

where = 3\/3ag(t^ - tl)/2. Since we are interested 
in nniaxially strained graphene along the zig-zag direc- 
tion in which ti = ^3, consequently c = 0. In this case. 
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FIG. 3: (Color online) Dirac point shift in respect to the 
noninteracting case as a function of the electron density ( 
in units of 10^^ cm~^) for strained and interacting graphene 
sheet. The Dirac point moves noticeably at the low electron 
density and tends to zero linearly in the high electron density. 



B. Interacting electrons in Hartree-Foclc 
approximation 

We consider the interaction of quasiparticles by using 
the leading diagram approximation which is the exchange 
interaction. In this sense, when interactions are treated 
in a mean-field approximation the Hamiltonian is written 
as H = i7o + Hmf whereP 

Hmf = - ^k-k'PcpC^'HLtPk^ , (14) 

kk' al3 

factor 2 stands for spin degeneracy and p is a den- 
sity operator given by Pai3{k) = ("-^°^l'0d">(V^([l + 

n?^ (■'/'o' l)ct/3 in which ni*^ are noninteracting band 
occupation factors, Vq is the Fourier transformation of the 
interparticle interaction. \tpo), the eigenvectors of non- 
interacting Hamiltonian when the full /(fc) expression is 
considered, are given by 



vl = 3al{4tl - tl)/Ah^ and vl = ^altl/Ah^. From this 
equation, it is obvious that the cross section of the en- 
ergy dispersion in strained graphene is an elliptic shape. 
Accordingly, the electron velocity is direction dependent 
and there is an anisotropy in Fermi velocitj^^. Notice 
that the energy at (<;f^,0) and (0,(7f„) is equal to the 
Fermi energy, ep with fcp = ^Jim. Furthermore the elec- 
tron density reads as wk = Qf^Qf and we thus have 



(12) 



ep = hy/Vj;Vyy/mr 



where gp^ = fcpJ, qf^ = kFS~^ and S — y^Vy/vx- 

Note that the electron density will change when we 
apply strain to a honeycomb lattice because the primitive 
cell area changes due to the strain as il' = ri(l -|- e^x + 
Eyy) = n{l + e(l — v)) so that 



no 



(13) 



where Uq is the electron density in unstrained graphene 



2 \±e-"''f 



(15) 



where exp[i0/(fc)] = f{k)/\f{k)\. After substituting 
Eq. [15] into the pseudospin density matrix, it can be de- 
composed into the charge and pseudospin-density contri- 
butions 



Pap{k) = f+{k)Sap + f-{k)nf.dap 



(16) 



where the unit vector is hf = x cos{(j) f) ~ y sm((j) f) and 
we have introduced a short-hand notation f± = {n'fj^ ± 

"f)/2. 

It is easy to see that the Hamiltonian consists of 
a momentum-dependent pseudospin effective magnetic 
field which acts in the direction of momentum k. The 
band eigenstates in the positive and negative energy 
bands have their pseudospins either aligned or opposed 
to the direction of momentum. Therefore, the interaction 
mean-field Hamiltonian changes to 



HhF = ^Ipiai SapBoik) + aai3.B{k)}lpkl3 

kaji 

J^-k-k'Uik') 

/flP ^ ^ ^ 

—vj:_g{xcos{ct)f{k')) - ysin(0/(fc')))/-(fc') (17) 

I 

where the integral term of _B(fc) is the exchange field and the sum of the first two terms are the band-structure 
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pseudospin magnetic field. It should be noticed that 
Eq. (17) is our main theoretical result. At zero tem- 



perature and for an n-doped graphene layer, f+{k) — 

^9{€f — Cfc) + ^ and /-(fc) — —^di^j: — (-f) where ep and 
eg are the Fermi energy and the noninteracting energy 
dispersion, respectively. To account partially for screen- 
ing and to avoid the well-known Fermi velocity artifact 
of mean-field theory in systems with long-range inter- 
actions, we have used an interaction potential including 
Thomas-Fermi screening 



Vk-k' 



e{q^^ + \k-k'\) 



(18) 



where q^^ = oteekp is the Thomas-Fermi screening vector 
and ckee — e^/iehvp) is the graphene coupling constant, 
which is density independent, e is an average dielectric 
constant of the surrounding medium. 

Notice that in our calculations, all integrals in the 
phase space are taken over the first Brillouin zone that 
is, a primitive unit cell of the reciprocal lattice ( see 
Fig. [l] (bottom). For the sake of simplicity, it is con- 
venient to work in a frame (^1,^2) where ki and ^2 
are along b'l and —6' 2, respectively. Therefore, it is 
accomplished by changing k^ — ^{ki + fc2)/|6'i| and 
ky = hi (ki — k2)/\b'i\ where the Jacobian is 



J 



6V3 



(l-^£)(l + £) 

9(1 - i^e)2 + 3(1 + e)2 



(19) 



The real and imaginary part of /(fc) are given by 
^ef{k) = — J2i ti cos(fc-(5;) and 9m/(fc) = — ti sin(fc- 
61), respectively. After having all ingredients, we are able 
to write the energy dispersion in HF approximation 



E±{k) = Bo(fc) ±^jB,ik)^ + By{k) 



(20) 



Once BQ{k) and B{k) are obtained, the ground state 
quantities can be calculated. Due to the exchange in- 
teraction, the Dirac points are shifted and a new Dirac 
point, Kjj is found by searching the zeroth of Eq. (20 1. 



C. Anisotropy renormalized Fermi velocity and 
charge compressibility 

The effective Fermi velocity is an important concept 
in Landau's Fermi liquid theory since it provides a di- 
rect measure of the many-body interactions in the elec- 
tron system. The low-energy electronic excitations in 
graphene are described by a massless Dirac Hamiltonian 
that is able to explain many transport properties. The 
Fermi velocity is the only parameter that appears in the 
model Hamiltonian and it plays the same role as the ef- 
fective masJ2ilin the standard Landau's Fermi liquid the- 
ory. The renormalized Fermi velocity can be calculated 
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FIG. 4: (Color online) Strain dependence of the renormalized 
velocity at (qf^ , 0) point, scaled by that of a noninteracting 
velocity as a function of the electron density ( in units of 10^^ 
cm~^) for different strain values and Oee ~ 0.9. Since the 
dispersion relation has an elliptic shape, the Fermi velocity 
has also an anisotropic form. For the positive strain val- 
ues, stretching the sample, the renormalized Fermi velocity 
at (qp^ , 0) point suppresses however, the renormalized Fermi 
velocity enhances by compressing the sheet. Importantly, the 
strain dependence of the Fermi velocity behaves differently as 
the sign of the strain is changed. 



when the ground state energy dispersion as a function of 
k is obtained. The renormalized Fermi velocity can be 
expressecP^ in terms of the wavevector derivative of the 
charge carrier dispersion relation evaluated at the Fermi 
surface 
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:^m\k=K 



k=Ko+qF 



(21) 



and therefore, the Fermi velocity component along zth- 
dircction is given by 



dB^ik) B^ik)'^ + By{k) 



hv* ^ 



dki 



dk, 



^B^{kY + By{kf 



\k=Ko+qF 



Once the ground state is obtained the compressibility, 
K can be easily calculated from 



1 



dn 



(23) 



where ^ is the chemical potential of quasiparticles and 
incorporates the kinetic energy and the first order ex- 
change interaction in the presence of the uniaxial strain. 
Therefore, it is easy to get the following expression. 



dn 



dk-p 

dn dk-p 



(24) 



in which dkp/dn = kp/{2n). Notice that the compress- 
ibility of the noninteracting system at low- energy is 
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FIG. 5: (Color online) Strain dependence of the renormalized 
velocity at {0,qFy) point, scaled by that of a noninteracting 
velocity as a function of the electron density ( in units of 
10^^ cm~'^) for different strain values and aee ~ 0.9. For the 
positive strain values, stretching the sample, the renormal- 
ized Fermi velocity at (0, qf^) point is enhanced however, the 
renormalized Fermi velocity is suppressed by compressing the 
sheet. The strain dependence of the Fermi velocity behaves 
differently with those results obtained at (gF^,0). 

given by Ko = 2/ (nh^VxVyy/Tm) . The expression reduces 
to the well-known expression of the noninteracting com- 
pressibility by setting = Vy. 




III. NUMERICAL RESULTS 

In this section, we present our calculations for the 
ground-state properties of graphene in the presence of 
uniaxial strain which we model as mentioned above. The 
shift of Dirac points, the anisotropy of the Fermi veloc- 
ity renormalization and inverse compressibility l/(n^K) 
are calculated by using the theoretical models and we 
compared the latter with the recent experimental mea- 
surements. In all numerical calculations presented here, 
we take into account the full Brillouin zone by using the 
expression of f{k). In this case, accurate calculations re- 
quire dense k— point sampling near the Dirac point and 
thus we use a dense adaptive sampling model the same 
as used in Ref. [33]. We have also dropped the constant 
1/2 contribution from the integrand of Bo{k) and the 
reason is that the latter term produces a constant energy 
contribution. 

As we have discussed in the previous section, the Dirac 
point moves to a new position by taking into account 
the full f{k) in strained and interacting graphene. The 
contour plot of the massless Dirac excitation energy is 
shown in Fig. 2 for aee = 0.5. Notice that the center 
of the Fig. 2(b) indicates the position of K^, ( the Dirac 
shift of strained graphene without the electron-electron 
interaction) and clearly Kjj moves to the left. It is also 
clear from the figure that the Dirac points, defined by 



FIG. 6: (Color online) The inverse compressibility [n^ft]"^ = 
dfi/dn (in units of meV 10~^°cm^) as a function of the elec- 
tron density (in units of 10^^ cm~'^) for different strain values 
for (a) a noninteracting and (b) an interacting massless Dirac 
electron system in graphene sheets. The symbols are the ex- 
perimental data by Martin et alW^. The strain dependence 
of the charge compressibility behaves differently as the sign 
of the uniaxial strain is changed. The coupling constant is 
aee = 0.26. 



E^{Kjj) = E_{Ku), are not exactly on the if^ and 
points. It is easy to show that By vanishes along x direc- 
tion for the uniaxial strain. Since Bo{K]j) does not van- 
ish, we might define an energy shift so that to get vanish- 
ing dispersion relation at the Dirac point. Therefore, the 
excitation energy is redefined as £{k) — E{k) — E{Kd). 
A deviation from the linear dispersion relation associated 
to the trigonal warping effect is clearly shown in Fig. 2(a) 
using the full expression of ,f{k). 

It is noticed that in a noninteracting graphene system, 
the Dirac point is given by Eq. ^ and it differs with 
K — 47r/(3\/3ao) when the lattice is unperturbed. Due 
to the charge-charge interaction, the Dirac point changes 
and the shift is found by searching the zeroth of Eq. ( 20 1 . 
The Dirac shift, Ko — scaled by k-p, as a function of 



the electron density ( in units of 10 cm ) is depicted in 
Fig. [3j The shift is noticeable at the low electron density 
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region and tends to zero linearly in the high densities. 
It is worthwhile mentioning that the Dirac shift value 
is enhanced when a graphene is compressed however, it 
is suppressed when a graphene sheet is stretched. In 
addition, the shift disappear^'^" in the linear limit of f{k) 
expression. 

Our results for the strain dependence of the Dirac elec- 
tron velocity, v*/vf, at {qF^,0) point as a function of the 
electron density n are demonstrated in Fig. 4. Since the 
dispersion relation has an elliptic shape, the Fermi veloc- 
ity has an anisotropic form. The numerical results are 
shown for different strain values, e and for the graphene 
coupling constant, age — 0.9. For the positive strain val- 
ues, by stretching the sample, the renormalized Fermi ve- 
locity at (qf^ , 0) point suppresses, however the renormal- 
ized Fermi velocity enhances by compressing the sheet. 
These behaviors are based on the effective hopping terms 
in the presence of the strain. Besides, the Fermi velocity 
decreases by increasing the electron density due to the 
exchange interactions between electrons near the Fermi 
surface. The strain dependence of the Dirac electron ve- 
locity at (0,(7Fj,) point, on the other hand, behaves dif- 
ferently with those results obtained at (gp^ , 0) point and 
our numerical results are shown in Fig. 5. As the uni- 
axial strain increases along the dirction of graphene de- 
formation, the Fermi velocity in this direction decreases 
quickly, whereas the Fermi velocity perpendicular to it 
increases. 

Recently, the uniaxial deformation of graphene on the 
unidirectionally modulated SiC steps shaped substrate is 
studied by using angular resolved photoemission spec- 
troscopy's!. A reduction of the Fermi velocity along 
the direction of graphene deformation, was observed, 
while in another direction the changes were not signifi- 
cant. Our numerical results on the Fermi velocity along 
the direction of graphene deformation are in good agree- 
ment with the recent observations. The reason for the 
insensitivity of the Fermi velocity perpendicular to the 
graphene deformation, Vy is probably due to the fact that 
the sample is shaped as a supperlattice along the direc- 
tion of deformation. The sample can be roughly realized 
as a combination of a pristine supperlattice graphene to- 
gether with a uniaxial strained ribbon. For the pristine 
supperlattice, Vy decreases and remains unchanged. 
On the other hand, we show that Vy increases when 
graphene is stretched and therefore there should be a 
cancellation between two contributions and it leads to 
insensitivity of the Fermi velocity perpendicular to the 
direction of graphene deformation in agreement with the 
experimental measurements. 

The local compressibility of graphene, on the other 
hand, has been measured "^^ using a scannable single elec- 
tron transistor. From the theoretical point of view, the 
compressibility was also studied by different groups at 
zero temperatur j^^l^^* ^ and also at finite temperatur^^. 
In order to calculate the charge compressibility in the 
presence of strain, we do need to calculate the deriva- 
tive of the total energy as a function of the electron den- 



sity. It should be noted that dBQ{k) / dk-p and dB{k) / dk-p 
change slightly by changing (j) — ta.n^^ (ky/kx) which is 
a consequence of the elliptic shape of the Fermi sur- 
face and the dependence tends to zero for small e 
values. Our numerical compressibility is obtained by 
taking an average over (j) and the results are shown 
in Fig. [6j The strain dependence of the compressibil- 
ity for a noninteracting system is also demonstrated in 
Fig. [6]ja). For the noninteracting case, it is easy to find 
that y/VxVy oc ?;f(1 — 9.1(1 — i^)e/(8to)) and thus for e > 
it gives rise to a reduction of the inverse of the compress- 
ibility. Importantly, the charge compressibility decreases 
(increases) by compressing ( stretching) a graphene sheet. 
Also the asymmetric response to the positive and nega- 
tive strain is evident in the charge compressibility as is 
evident in the renormalized Fermi velocity. 

In Fig. we compare our theoretical predictions for 
the inverse compressibility of doped graphene with the 
experimental results of Martin et aZ.'^-'as a function of the 
electron density in units of 10^^ cm~^. Martin et aL^SS fit- 
ted the experimental inverse compressibility, {n?K)~^ to 
the kinetic term using a single parameter Fermi velocity 
which is larger than the bare Fermi velocity. Note that 
the kinetic term in graphene has the same density depen- 
dence as the leading exchange and correlations termJl^l 

As it is clear in Fig. |6][b) the inverse compressibility 
of an interacting system is higher than the experimental 
value. By increasing the uniaxial strain effects, i.e., in- 
creasing the constant strength, e our theoretical results 
move down. Therefore, including the exchange and strain 
effects in our theory, gives results very close to experi- 
mental data. 



IV. CONCLUSION 

We have studied the ground state properties of a 
graphene sheet within the Hartree-Fock theory incorpo- 
rating the uniaxial strain in the system. We have shown 
that the quasiparticle properties could be strongly strain 
dependent and substantially different than the usual pris- 
tine graphene sheet. The Dirac points move due to the 
strain and the anisotropy of Fermi velocity renormal- 
ization is obtained. The Dirac electron Fermi velocity 
is highly dependent on the sign of strain ( stretching 
or compressing) even if the strain is negligibly small. 
The renormalized electron velocity along the direction 
of graphene deformation decreases with stretching how- 
ever, it increases by compressing the graphene sample. 
The Fermi velocity along a direction perpendicular to 
graphene deformation behaves in a different way. We 
find the reduction of Fermi velocity along the direction 
of graphene deformation to be in good agreement with 
the recent experimental observation. 

Our calculations of inverse compressibility compared 
with recent experimental results of Martin et alW^ 
demonstrate the important influence of strain on the 
thermodynamic quantities in a graphene sample. 
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We remark that in a very small density region, the 
system is highly inhomogeneous and the effect of strain 
might be very essential. A model going beyond the 
Hartree-Fock approximation is necessary to account for 
increasing correlation effects at low density. 
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